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^ . Abstract 
Q 

. Spiral waves arise in many biological, chemical, and physiological sys- 

tems. The kinematical model can be used to describe the motion of the 
spiral arms approximated as curves in the plane. For this model, there ap- 
peared some results in the literature. However, these results all are based 
upon some simplification on the model or prior phenomenological assump- 
tions on the solutions. In this paper, we use really full kinematic model to 
classify a generic kind of steadily rotating spiral waves, i.e., with positive 
^ ' (or negative) curvature. In fact, using our results (Theorem |5J, we can 

answer the following questions: Is there any steadily rotating spiral wave 
for a given weakly excitable medium? If yes, what kind of information 
we can know about these spiral waves? e.g., the tip's curvature, the tip's 
tangential velocity, and the rotating frequency. Comparing our results 
with previous ones in the literature, there are some differences between 
them. There are only solutions with monotonous curvatures via simplified 
model but full model admits solutions with any given oscillating number 
' of the curvatures. 
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1 Introduction 

Rotating spiral waves appear in many biological, physiological, and chemical 
systems. For example, spiral waves arise in cardiac arrytmias |2(JI I3UI 13*21 1461 16 1| . 
egg fertilization the Belousov-Zhabotinsky (BZ) chemical reaction |4*l 1251 
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I5T1 13*71 IT71 """""1 177) . catalysis ["03 02], and aggregation of slime mold amoeba 
P E21 ESI EH]- These waves have been studied extensively from experimental 
jU |2Sl E3 S3 IIHI E21 ESI EU EZI, numerical [231 [2^1 [271 HM FTol 17^1 [75]. and 
analytical H21 EH E21 EH EHl EBl E21 E31 EH EH1 EH] aspects. 

Often, such systems exhibit the so-called "excitable" property A 
spiral wave represents an excited moving spatial region with two spiral-like, thin 
phase-transition layers (or interfaces): wave front and wave back. These layers 
separate two different phase regions, excited and refractory (or recovery), and 
intersect at a tip (see Figure 12.5 in [""""j or Figures 13, 17, and 22 in 40 ). 

Usually, the appearance of spiral waves is not desirable. For instance, heart 
attacks are caused by an abrupt shift from rhythmic pumping to spasmodic 
convulsion of the heart. Normally, with each heartbeat, an electric wavefront 
propagates across the interconnected muscle fibers, causing them to contract. 
However, because of some abnormality in the tissue, this wave can become stuck 
and start rotating as a spiral wave [SHUSH]- Furthermore, if spiral waves occur 
in the cortex, they may lead to epileptic seizures. On the retina or the visual 
cortex, they may cause hallucinations [8] E] • 

Spiral waves can be created by breaking propagating waves [23 EDI ESI- I n 
|60| . a circular wave is broken mechanically by ejecting a gentle blast of air 
into a small section of the wave. The process by which a broken wave evolves 
toward a rotating spiral wave has been addressed in []"""] [""H] EH El EI EE] . As 
pointed out in |19| , along the interface separating the excited and the refractory 
regions, there exists a point where the normal velocity of the interface changes 
sign. This induces a twisting action on the interface motion, and as a result, 
the interface starts wrapping around some center of rotation, to form a spiral 
structure. Eventually, in an appropriate parameter range, a state of steady 
rotation is achieved, during which the spiral tip traces a circular trajectory 
at constant angular velocity, and the curvatures of the spiral-shaped interfaces 
keep constant in time. Such steadily rotating spiral waves have been observed 
in many experiments [481 1491 ITT) . 

Another way to initiate spiral waves is by applying a spatially graded per- 
turbation to a medium in the recovery phase after excitation [721 173| . Besides 
steadily rotating spiral waves, non-steady forms of rotation, known as "mean- 
der" US], are also confirmed by experiments [251 1601 IB!"] and numerical simu- 
lations [S] ESI ESI EH EE] ■ When the controlled parameters are changed within 
some suitable range [^SHSSIES], or due to the anisotropic influences of the media 
|61| . or the effect of electric field ^50, EI], a transition from the steadily rotating 
spiral waves to meandering ones may occur. To study the above mentioned 
wave patterns, many models have been proposed, which we can categorize here 
as "phase field" type and "sharp interface" type models 1 . 

For the phase field type of models, well-known examples are the reaction- 
diffusion equations (for example, the FitzHugh-Nagumo model for the action 
potential wave in neurons [SO]), the Oregonator model for the BZ reaction [""|"J 
16 7| . the bidomain model |SS], and the chemotactic models (SSJESI- The bidomain 

x We borrow these two terms from material science. 
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model is a mixture of elliptic and reaction-diffusion type of equations. It is used 
to describe the electrical activities of cardiac tissue. The chemotactic models 
describe the oriented movement of cells in response to the concentration gradient 
of chemical substances in their environment, such as the aggregation phenomena 
of slime mold amoeba. 

Using the reaction-diffusion equation, the existence of steadily rotating spiral 
waves near a homogeneous steady state has been proved formally by Hagan |^] 
and later rigorously by Scheel [HI]. Other more complicated wave patterns, e.g., 
meandering and hypermeander were also explored by El E21 EHl HH ■ 

For the sharp interface type of models, some typical examples are the free 
boundary models Q71 H3 El EH EID 023 ESI and the kinematic models [H 
HS1 E3 El US EH Other models are given in [HEH7S] . 

If the two phase transition layers, as mentioned above, are very thin, we can 
view them as two plane curves meeting at the tip. Then, the free boundary 
model consists of an interface equation, a simplified phase field equation, and a 
dispersion relation with some boundary conditions. The interface equation (13) 
or (38) in jSS] is derived from the so-called eikonal-curvature equation (11) there 
and governs the motion of a wave front; the simplified phase field equation (37) 
in [HE] is derived from the original phase field model and affects the velocity 
appearing in the interface equation; and the dispersion relation (see section 3.3 
in [HHj) reflects the interaction between spiral waves (refractory tail effect). 

In weakly excitable media, spiral waves rotate around a large circle, the 
excited region becomes rather narrow, and the spiral waves are sparse. In this 
case, we need to consider the wave front only. Thus, the spiral wave can be 
viewed as a single curve, which performs a motion in the plane regardless of the 
interaction between the spiral waves. That is, the dispersion relation can be 
ignored. 

The kinematic model is then formulated in terms of a motion of a single curve 
with free end; see section 2 for the details. The validity of this kinematic model 
has been verified, since various wave patterns appearing in the experiments can 
also be produced with this approach [44] . 

Finally, we should point out that many main governing equations of the 
sharp interface models can be derived from the phase field models by taking 
various types of singular limits, e.g., the small parameters may represent the 
width of the thin transition layers. Therefore, it is natural to expect that there 
should be a close relation between the solutions of these two type of models. 
Indeed, it has been found that there is a good agreement for temporal periods 
of rotating spiral waves, when comparing the reaction-diffusion approach with 
the kinematic one |45|. 

In this paper, we use the really "full" kinematic model to prove the existence 
of steadily rotating spiral waves. In fact, using our results (Theorem[HJl, we can 
answer the following questions: Is there any steadily rotating spiral wave for a 
given weakly excitable medium (i.e., Vq, D are given below). If yes, what kind of 
informations we can know about these spiral waves? e.g., the tip's curvature («o 
below), the tip's tangential velocity (G below), and the rotating frequency (u> 
below). In the literature, there appeared some results about kinematic model. 
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However, these results all are based upon some simplifications on the model 
or prior assumptions on the solutions. For examples, [T2|, E2> an d F^l 
considered only simplified model (see the Remarks in Section 2); |23], |44|. 
and [ZH] made a prior assumption on the solution, i.e., tip's tangential velocity G 
is 0, which was proven to be a very special case of steadily rotating spiral waves in 
|T3) : pi) is valid only for nonrotating ones, i.e., angular frequency ui is (see the 
explanation after lfTB|l below); |2U considered solutions with positive curvature 
k. In this paper, we also consider solutions with unchangeable curvature's sign 
(i.e., with positive or negative curvature). However, this condition is generic. 
Note that the formula (3.11) in [23 is w(0;0, 1(0)) in Theorem|Hl i.e., a special 
case of our results with zero tip's tangential velocity (G = 0) and monotonous 
curvature (i — 0). Note also that steadily rotating spiral waves with positive 
curvatures must have monotonous curvatures via simplified model (see \'21\): 
however, using full model in this paper we have confirmed the existence of 
oscillating ones, i.e., the case "i ^ 0" in Theorem|Hl 

In section 2, we recall the kinematic model following Mikhailov et al. 02 
and a new version by ^H]. In section 3, we state and then prove the main 
theorem of this paper. Numerical results are presented in section 4. 

2 The Kinematic Model 

The kinematic model, formulated in terms of a motion of curves, was first pro- 
posed for describing excitation waves in a cardiac muscle by Wiener and Rosen- 
blueth |7()j . They assumed that a plane curve moves in a normal direction with a 
constant velocity. They showed that such a curve, rotating around an obstacle, 
forms a spiral (which represents an involute of this obstacle) and approaches an 
Archimedian spiral far from it. 

Since there occurs a singularity in the curvature of this spiral curve near the 
tip jZO]) and also by the subsequent analysis of wave propagation based on an 
excitable reaction-diffusion equation, it has been found that the normal velocity 
should be relevant to the local curvature [23 OH EEJ- For example, the so- 
called eikonal- curvature relation or mean- curvature flow type equation (see 
below) has been shown to be correct when linear approximation is considered 

[UlEllEH]. 

We should note here that Wiener and Rosenblueth considered only pinned 
spirals that rotate around an obstacle. However, subsequent numerical simula- 
tions and experiments with various excitable media have revealed that the same 
media can also support spiral waves free of any obstacles |fc>51 172] . 

Later, Mikhailov et ai.|33] proposed the following standard kinematic model: 
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dn Ok ( f s \ d 2 V 

~di + lh[ KVd t + G ) +k2v + -Qgl = °> 0<s<oo, 0<t<oo, (1) 

V = V - Dk, < s < oo, < t < oo, (2) 

G = Go-rn , < i < oo, (3) 



dK _ q® k 
dt ds 



, < t < oo, and (4) 

s=0 

lim k(s, t) = 0, < t < oo, (5) 



where k = k(s, t) is the curvature of a plane curve depending on time t and 
the arc length s measured from the free tip, V = V(s, t) is the normal velocity, 
G = G(t) is the tangential velocity of the tip, kq = :— lim s ^o n{s,t), and 

Vo, D, Go, r are parameters determined from the media. 

Here, QJ is a general equation satisfied by a plane curve of any point at 
which the curve moves in the normal direction (see (4.41) in 0Oj or (2.7) in 
|44p; (0 is the eikonal-curvature relation, which relates the normal velocity of 
the curve to the local curvature (see (4.10) in 0U| or (1.1) in 44 ); © is the 
tangential velocity of the tip 2 and |0J and JSJ) are boundary conditions (see 
(2.10) in gl|). 

We explain J2J) as follows. The curvature near the tip, hence Ko, affects the 
normal velocity near the tip by @. Then, by (3.27) from the normal veloc- 
ity near the tip affects the width of the excited region near the tip. Moreover, 
it has been observed that the width of the excited region near the tip influences 
the tangential velocity of the tip (23 0^. Hence kq is related to G. 

Notice that (j2J and (JHJ are good linear approximations of the perturbation 
expansions only when the curvature is small. Hence "k <C W in © (or v ~ vq) 
is required. 



Remarks Elkin et al.\T$\ used perturbations techniques and Fredholm alter- 
native theorem to obtain more general equations compared to © and Q) here 
(see (28), (32), (36), and (37) in [TijV In fact, it has been found that © is 
only a special case of (32) in ^2] (see the context below (37) in ^1]) and © 
is not correct. Later Elkin et al. ^H] solved the above new model under some 
simplification, i.e., replacing V with Vb in the term "J" s nVdl; " in © (see (4) 
in [SI). 

In this paper, we consider steadily rotating spiral waves (so §f = 0). Let 
G be a constant and let co be the momentary rotational angular velocity of the 
tip. Then, using J2J and integrating the both sides of © give rise to: 

«(*){ f KOiVo - Dk(OH + G}- Dk'(s) = u. (6) 
Jo 

2 12J is postulated in |B]|7| and derived by perturbation techniques in Appendix A of 1441 . 
when linear approximation is considered. 
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Thus it follows that 

«'(0) = ~*(0) - g = - (7) 

(see also (6) in JE])- After differentiating © once and replacing the integral 
term, we obtain: 

-Dk" + k 2 (V - Dk) + — (lu + Dk') = 0, (8) 

provided k^O. Without loss of generality, we assume k > hereafter, since if 
n < 0, then Q and (jSJ) are invariant by replacing k, w, Vq with —k, —u>, —Vq, 
respectively. Let l(s) := ln«;(s). Then JSJ), Q, and © give rise to: 

^._ 9(i(s)) + (9) 

i(0) = lnK , (10) 

<'<»' = -7^ + §• 

lim 1(a) = -oo, (12) 

s — >oo 

where <?(Z) := e 2 ' — e z . The above equations @, (fTTTf) . (fTTf) are also equivalent 
to the following ones: 

l'(s)=v(s) (13) 
v'(s) = -g(l(s)) + ^e^v(s) (14) 

u( 0) = -£e-<(°) + g. (15) 

Furthermore, if w(s) ^ 0, then the solution v(s) in (|13|) . l|14l) can be viewed as a 
function of I by the inverse function theorem and satisfies the following equation: 

* £<& + £e-'. (16) 
at v D 

Note that in |23], Ishimura et al. considered the above problem with 
G = and l'(0) = (see (13) in 24 ). Therefore, the results in are valid 
only for the case "u> — 0" . 

Now for any excitable media with given parameters Vq > 0, D > which are 
fixed throughout this paper, we want to find suitable G, to, and ^(0) (i.e., /to) 
such that Q12JI. H13(l . 1)14(1. and (|15|l has global existence solutions. Let £?(Z, u) := 
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3 Behavior of the Solutions 

Lemma 1 For any to < 0, there are no global existence solutions of ^13\) , \Tljty 

satisfying lim l(s) — — oo. 

s >oo 

Proof. Since lo < 0, given any solution of (fb^ . fTljl . say, (l(s;ui), v(s;uj)) (drop- 
ping out the dependence on initial data here for simplicity), E(l(-;u)),v(-;w)) 
is a decreasing function by (11) of [1]. Note also that E(l,v) > S(ln-^,0) = 
^r(-rr) 2 - Thus E(l(-; lo), v(-} lo)) is bounded. If l(-; to) exists globally with lim l(s 

— oo, then we have lim se 21 ^ = ^- by LEMMA 3 of H3|. This is impossible, 



.4 — -oo 



since lim se 21 ^ > > t£t. This completes the proof. ■ 

s-»oo — 2V 

By above lemma and THEOREM 2 of we only need to consider lo > 



Lemma 2 Given any global existence solution of J ^T^\j , (l(s;u)),v(s;u>)), if 
v(0;lo) > 0, then there is some a > swc/i £/ia£ v(a;Lo) = (no matter the sign 
ofui). 

Proof. Suppose not, i.e., (l(s;u>),v(s;u>)) always stays in the upper half plane: 
v > 0. Note that the upper half plane is splittcd into three parts: R\ := {(I, v) : 
v > and -g{l) + %e~ l v > 0}, R 2 := {(l,v) : v > and -£f(Z) + %e~ l v = 0}, 
and i? 3 := : « > and -g(l) + ^e~ l v < 0}. If (1(s;lo),v(s;lo)) always 

stays in R%, then both l(-;u>),v(-;u) increase as s — > oo by l|13fl . (|14fl an d so 
both lim 1(s;lo) and /x := lim v(s;lo) > exist. Thus clearly, we have lim 

s — >oc s — >oo s — >oo 

l(s; lo) = oo. Since v (s; to) > and lim —g(l(s; lu)) — — oo now, we have \i = oo 

s — >oo 

by (|T1J| . On the other hand, it follows from ((Til) that: 

«(»; lo) + ^ e - 1 ^ ~ (v(0; lo) + ^e~ 1 ^) = - J* g(% Lo))dt. (17) 
Taking limit as s — > oo on the both side of (|17f) . we have 

/•OO 

lim [„(«; w) + ^e"^) - (v(0; lo) + ^e~ 1 ^)} = oo = - / ff (Z(f ; W ))dt. 

(18) 

However, — g(l(t;Lo))dt = — f Q g(l(t;co))dt — g(l(t;Lo))dt, where 1(S;lo) — 
ln-j^. Since g(Z) > as 1 > In and l(t;co) f oo as t / oo, we have 
— <?(Z(t; w))c?t < and so 

g(l{t:Lo))dt < g(l(t; u))dt < oo. (19) 

j o 

Thus (|19(l contradicts (|18|l . Now the only possible case is that (1(s;lo),v(s;lo)) 
goes through the curve i?2, enters the region i?3, and then stays there for ever. 
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Hence 1(-;oj) and v(-\uj) are increasing and decreasing functions respectively 
now by lfT3 Tl .lfT3 1l . Clearly '■= lim v(s:lj) 6 [0,oo), since v(s;lj) > 0. 

s — *oo 

Note that the "Poincare-Bcndixson theorem" implies lim 1(s;lj) — oo; oth- 

s — >oo 

erwise, (Us; uj), v(s; u>)) will stay in a bounded region and so lim 1(s;lo) — 

s — >oo 

In and lim v(s\lj) — 0, since (In -yii 0) is the only possible lj— limit point of 
(l(s; lj), v(s; lj)). A contradiction then occurs, since lim 1(s;uj) ^ In ^ by (HHJ) . 
Now it follows from " lim l(s; lj) — oo" and (|14|l that lim v (s; w) = — oo which 

s — >oc s — >oc 

contradicts "v(s;lj) > 0,Vs > 0". ■ 



Lemma 3 For lj > 0, given any solution of \1!^) , ^T^ , (1(s;lj),v(s;lj)), with 
maximal existence interval [0,r), if E(l(s^; lj), v(s#; lj)) > for some S 

[0,r), then r < oo and lim 1(s;lj) = lim v(s;lj) — — oo, i.e., blows up in 

s — >r~ s — 

finite s. 

Proof. Assume the contrary that r = oo. Then by Lemma|21 (1(s;lj),v(s;lj)) 
will enter the lower half plane after a short time, say, a. Since E(1(-;lj),v(-;lj)) 
is increasing and E(1(s^;lj),v(s : ^;lj)) > for some S [0,T), we have 
E(1(s;lj),v(s;u))) > Vs e (s#,oo). It then follows that v(s;lj) < 0, Vs e 
(max{cr, s^}, oo), since the level curve E(l, v) = is a barrier preventing (l(s; lj), v(s; lj)) 
from touching the I— axis. Clearly, both \x\ := lim/(s;w) < oo and /12 '■= 

s — »oo 

lim v(s; lj) < exists. Note that [ii ^ 0; otherwise, again by the "Poincare- 

s — >oo 

Bendixson theorem" we have fi\ = — oo and so lim E(l(s; lj), v(s; lj)) = 

s — >oc 

which is impossible. Now it follows from "/12 G [— 00, 0)" that [i\ = —00. By 
1)14(1 . we then have lim v {s;u>) = —00 and so /i2 = —00. Moreover, we also have 

s — >oo 

lim v (s;lj) = —00 by differentiating both sides of (|14|l and taking limit. 

s — >oo 

Multiplying both sides of H17|) by e 1 ^'^ and integrating them, we then have 

£(s;a>) 



^+/ e l dl-(v{0;Lj) + ^-e- l(o ^) I e l ^dp = -l e l ^l g(l(t;uj))dtdp. 



(20) 

Now we want to show that e 1 ^'^ = o(^) as p — > 00. Utilizing PHopital's rule, 
we have 

lim ^— = lim — r, r = lim ■ ■ r , 7 = lim 



p^oo p 3 p— too e~^ p ' w ' p^oo — v(p; w)e _ ^ p; ") p^co — y' (p; Lj)e~ lt ^ p '^ + v 2 (p; Lj)e~ lt ^ p '^ 

(21) 

6 

— lim = 0. 

p^oo [— v"(p] uj) + 3v (p; u))v(p; to) — v 3 (p; co)]e~ l ^ p '^ 



Since 5 (0 > ff (In we have - / °° e l M J P w))«ftdp < -gQn ^) / °° e'to-Vp < 

rl(s;ui) 1 p — 00 



00 by (|2TJ- Also, lim Jhq^ e l dl = L^^dl E (-00,00) and J °° e l{p > u) dp £ 
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(0,oo) by l|21l) again. By taking limit as s — ► 00 on the both sides of (J2DJ, a 
contradictions occurs. Therefore T < 00 and so the proof is finished. ■ 



Lemma 4 Given any to S (0, -jy-), i/iere exists a unique I* = 1*(lo) <G (In In ^jf-) 
such that has a global existence solution (1(s;lo),v(s;lo)) defined on 

(—00, 00) with (1(0; lo),v(0; to)) — (I*, 0), l(s; lo) \ —00, v(s; lo) — * as s — ► 00, 
and v(s;lo) < for all s > 0. Moreover, (1(s;lo),v(s;lo)) rotates counterclock- 

Vo 
D 



wise around (ln^-,0) and tends to it as s decreases to —00 



Proof. By THEOREM 3, LEMMA 4, and LEMMA 5 in [1], we just know 
there exists an interval [sup A, inf B] C (In In =jy-] such that the solution 
of (|13() . 114l) with initial data (Iq,Q), where Iq € [sup A, inf B] , has the prop- 
erties mentioned in this lemma, denoted by (Z(s; u), Iq), 5(s; lo, lo))- Note that 
by Lemma we have [sup A, inf B] C (In -g-, In ^yy). Now we want to prove 
sup A = inf B. 

Since the solution v(s; u>, Iq) 7^ 0, for all s £ (0, 00), v(s; lo, Iq) can be viewed 
as a function of / by the inverse function theorem, i.e., 5(s; lo, Iq) = v(l; lo, Iq), for 
/ G [— 00, lo] and so satisfies Hlfi[) . For any two solutions v(1;lo,Iq\), v(1;lo,Iq2), 
denoted by v\(l) and 1*2 (0 respectively for simplicity, by l|16|l we have 

d ^^= g M{vi-v.). (22) 
dl V1V2 



Thus (5k - v 2 )(l) = (Si - 5 2 )(ln exp(- f/ n ^ o!Z) for / < In Sup- 



VI (1)V2(1) ' 

pose (v\ — 5 2 )(ln^-) ^= 0. Then since 5i,5 2 < and g(l) < for I < hiy|, 

we have lim (5i — 5a)(Z) 7^ O.This is impossible, since lim 5(s;w, £01) = 
; — >— 00 s — »oo 

lim v(s;lo,Iq2) = 0. By uniqueness, we have Iq\ — Iq2 and so sup A — infi?. 

s >oo 

Let 1*(lo) := supA = inf B and (1(s;lo),v(s;lo)) := (1(s;lo,1*(lo)),v(s;lo,1*(lo))). 

Finally, since the a— limit set of (1(s;lo),v(s;lo)) consists of only one point 
(ln^,0) which is a unstable focus, (1(s;lo),v(s;lo)) is also defined on (— 00, 0] 
and lim (I(s; w), 5(s; to)) = (In ^, 0). ■ 

s — > — 00 

Given any solution of lfH5|) . lfl4^l . say, (Z(s; co),v(s; lo)), E(1(-;lo),v(-; lo)) is an 
increasing function, since lo > 0. Also, by LEMMA 5 in [52], (1(s;lo),v(s;lo)) 
always stays inside the curve: E(l, v) = in the I — v phase plane. 

Proposition 5 Suppose lo £ (0, -jy-). for any global existence solution of ifi3j) , fm ), 

iis or&«i curve in the l — v plane coincides with that of (1(s;lo),v(s;lo)), i.e., they 
both lie on the same solution curve in the I — v plane. 

Proof. Let (l(s; lo), v(s; lo)) be a global existence solution of l|13[l . (|14() . i.e., existence 
interval contains [0,oo), (also dropping out the dependence on initial data here 
for simplicity). Suppose (1(s;lo),v(s;lo)) stays in a bounded region in the / — v 
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plane. Since the equilibrium (lna,0) is a unstable focus, the only possible 
lo— limit set of (l(s;co),v(s;ui)) consists of periodic orbits by the "Poincare- 
Bendixson theorem". However, by the "Bendixson criterion" (see [ ]), there 
is no any periodic solution of (|13J) . (|14|) . Thus our assumption is wrong. There- 
fore, any global existence solution is unbounded. 

Note that by Lemma[31 (i(s; to), v(s; lo)) always stays inside the curve:i?(7, v) = 
in the I — v phase plane. Since (Z(s; u>), v(s; lo)) is unbounded, we can choose 
some si large enough such that (l(s±\w),v(si\w)) is located in the lower "left 
tail" region : {(l,v) € R 2 \E(l,v) < 0, v < Q, and I < II, for some small enough 
with E(1(si;oj),v(si;oj)) > E(l*,0) by the structure of level curve of energy 
E{l,v) with negative value, i.e., E(l,v) — c < 0. Suppose v(s2',lo) = for some 
S2 > 8\. Then by the unboundedness of (1(s;co),v(s;uj)), the direction of flow 
on the upper half plane, and the increasingness of E(l(-;ui),v(-; ui)), there exists 
some S3 > S2 such that v(sa;uj) — and 1(s3;lo) G (l*,ln^ L ). Note that any 
solution curve which intersects the I— axis at (Ir,0) with Ir S (Z*,rn^-) must 
touch the level curve E(l,v) = and so blows up in finite time by Lemma 
since I* = inf B (see [21] for the definition of set B) and the uniqueness theo- 
rem of ordinary differential equation. Thus (1(s;lo),v(s;lo)) blows up in finite 
time which leads to a contradiction. Therefore our assumption is wrong. We 
conclude that v(s;lo) < for all s > si which leads to the following facts: 
lim s l(s; ui) = — oo and lim s too v(s; lo) = 0. 

Now using similar argument as in LemmaQJ the solution curve of (/(s; lo), v(s; lo 
must coincide with that of (l(s; lo), v(s; lo)) in the / — v plane. ■ 

Let 1* 1r (uj) := 1*(lo), 1* 2R (lo), l* R {uo), ... > In ^ and i* £ (w), 1* l (lo), ...,1* l (lo), . 
< ln^ be the sequence of intersection points of (1(s;lo),v(s;lo)) with the 
/—axis greater than and less than hi-^ respectively as s — > — oo. Note that 
1<k(<") > 1 * 3 r(u) and l* L (iu) < l* jL [io) as i < j. 

Lemma 6 (i). l* R {-) and l* L {-) are decreasing and increasing functions re- 
spectively on (0, —n~)> 

(ii). l* R (-) and l* L {-) are continuous functions on (0, — jy-). 
(Hi), lim 1* r {lo) = In ^jy- and lim 1* l (lo) = — oo, i = 1,2,3,... 

LO— >0+ Ul^0 + 

Proof. Since v(s;u>i) ^ Q,v{s;lo2) ^ 0, for all s € (0,oo), v(s; L0\), v{s; LO2) can 
be viewed as a function of I by the inverse function theorem, i.e., v(s;coi) — 
v(1;loi),v(s;lo2) = v(1;lo 2 ) for I E [— 00, 1* r (loi)), I € [— 00, II r (lo 2 )) respectively. 
Note that for lo\,lo2, the corresponding solutions Vi(l) := v(l;u>i), i>2(l) '■= 
v{l\L02) of l|16() have the following relation: 



(vi-v 2 )(l) = 
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K-, 2 )(r)exp( f -^-dt)+(^^)eM /' -#7^*) /'«p(-t- /' -M^dp)^ 
J T vi{t)v 2 {t) D J T Vi{t)v 2 {t) 7 r 7 r «i(p)u 2 (p) 

(23) 

whenever both v\,v 2 exist on [r, I] or [/,r]. Now apply (|2*3^l to «i(Z) := w(Z;wi), 
v 2 (l) :— v(l; UJ2), for Wi, w 2 6 (0, ^jy-) with a>i < cj 2 . 

First we consider ?J fl (-)< 

Case 1. /^(^i) < ^ fl (o> 2 ). 

Then there is a small £ > such that Vi(l),v 2 (l) are defined on (—00, l* R {<*Ji) — 
e] with Wi(Z*jj(wi) — e) > ^(^(cji) — e). Taking r = l* R (oJi) — £ and I < r, by 
(|23|) we obtain: 

(ui - u 2 )(0 > («i - «2)(t) exp( / L~ M dt). 

J T Vl{t)V 2 {t) 

Since lim («i - u 2 )(r) exp( f r v l (tyi 2 (t) dt ) > °> we nave lj m (£1 ~ ^ > °- 
A contradiction occurs. 

Case 2. l\ R (wi) = Z* i? (w 2 ) and there is e > such that v\(l\ R {u>i) — e) > 
#2(^(^1) - e). 

Then taking t = — e and / < r, by l|23|) we obtain: 

(-i-- 2 )(0>( i ^)exp(/ i _f^ 7 -^) /'exp(-t- fBtodfidt. 

D J T vi(t)v 2 (t) J T J T vi_{p)v 2 {p) 

Since ; UmJ^) exp(^ sj^dfc) £ exp(-t- J* ^_dp)dt > 0, we also 

have a contradiction. Hence the only possible case is ll R (u>i) — l* R (w 2 ) and 
Vi(l) < v 2 (l) for I G (—oo,II r (loi)) other than > l* 1R (uj 2 ). Taking any 

t G [ln-^,;*^!)) and J G [r, Z* H (wi)), by we have 

ivx-v 2 )(l) < (^j^)eM f ^^dt) feM-t- f ~Mr-^ dt < °- 

L> J r Ul(t)« 2 (t) J T Jr Vl(p)V2{P) 

(24) 

Clearly, the right hand side of l|24|) will not tend to 0, as I — > Z*^(k>i) . This 
contradicts to lim (vi — V 2 )(l) = 0, since 1* r (uji) = l* R (^2}- 

By above discussions, we conclude that l* R (u}i) > l* R (w 2 ), if < u 2 . 

Now we consider l* L (-)- Given any sl ± < such that the solution curves 
of (l(s;u>i),v(s;wi)) and (l(s; u 2 ), v(s; uj 2 )) stay in the upper half plane when 
s G [sij,0), both w(s;wi) and v(s;lu 2 ) then can also be viewed as functions of 
I, vi(-), v 2 (-) as before and so is also valid. Due to II r (uji) > 1\ r [lu 2 ), we 
can take r = ll R (ui 2 ) — 8 for some small 8 > such that both «i(-) and v 2 {-) 
are defined on (max{ZJ L (o; 1 ), ^ L (w 2 )}, r] with Vi(t) > v 2 (t). By we have 

(ft-W) > (^J^)eM f ^M^dt) feM-t- f ~MT-,d P )dt > 0, 

(25) 
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for I G (max{ij' t (wi), ll L (ui 2 )} 7 r). Since the right hand side of will not 
tend to as I — ► max{^ L (cji), /^(l^)} -1- , we conclude that < ^ L (w 2 ), 

if U>i < L0 2 - 

Arguing in the same way as above, we can derive l* R {tUi) > l* R (ui2) from 
"'(Vi)lM < l l-i)L^T and then l* L {u x ) < l* L {w 2 ) from «l* R {u x ) > l* R {u 2 f 
for i = 2,3,4, ... This completes (i). 

Now we want to prove (ii). First consider l* R (-)- Note that by (i) and Lemma 
|HZ := hm l* R (w) exists for any given uj q G (0, ^-) and l e (ln^,ln^]. 

UJ — >UJQ 

Suppose Iq = In i.e., the initial data and parameter pair {{l* R (u), 0), u>) 
tends to ((In 0), ujq). Then by the "continuous dependence on the initial 
data and parameter" theorem of ordinary differential equation and Lemma |3 
a contradiction occurs, since -E(ln 2gi,0) = and E(1{-;u),v{-;uj)) is increas- 
ing. Therefore we have Iq G (hi j), m ^d 1 )- Also by the "continuous depen- 
dence" property, the solution of H13l) . ill4[) with initial data and parameter pair 
((Z ,0),wo) must stay in the region {(Z,u) G R 2 \E(l,v) < and v < 0} . Then 
by LEMMA 1 in [1] this solution is global existence one. Thus by Lemma 0] 
and Proposition^! we have Iq = Z lfl (wo) and so lim II r (lu) = II r (luo). Hence 

U> — >LJQ 

/*(•) = l*[ R (-) is a continuous function on (0, -p 2 -)- 

For l* L (-),i = 1,2,3, ...and l* R {-),i = 2,3, by the continuity of ll R (-) and 
the "continuous dependence" property of o.d.e. when choosing long enough exis- 
tence intervals [s£ 2i ,0], [si 2i _!,0] such that the orbit curve of (Z(s; loq), v(s; ujq)) 
intersects the I— axis 2i and 2i — 1 times respectively, both 1* l {lo) and l* R (ui) 
will stay close to l* L (u>o) and 1* r (loq) respectively as oj is close to uq. Therefore 
^il(') anc ^ are continuous functions. 

Finally we prove (iii). Samely, the limit £ := lim 1* r (lu) exists with £ G 

(In §,ln Suppose C € (In ^, In Consider the solution of fE ty -fU Jl with 
initial value-parameter pair ((£, 0),0); it is periodic with some period, say, T G 
(0, oo). Applying the "continuous dependence" property of o.d.e. when choos- 
ing existence interval [0,iT], solution (1#(s;oj),v#(s;lu)) := (l(s + s i ( LU y,ui),v(s + 
s i(u>)] w )) over the existence interval [0, iT] where (Z(sj( w ) ; u>), w(sj( w ); u>)) = (1* r (lu) 7 0), 
must intersect the "negative" I— axis (less than In at least i times as ui — > + . 
This is a contradiction. Thus lim 1* r (uj) — m^. 

By (i) again, we have <r := lim l* L (w) exists and c G [— oo, In -^). Suppose ^ 

G (— oo, In ^-). Now we apply the "continuous dependence" theorem of o.d.e. to 
the solution of 1|13|) . 1)14(1 with the initial value and parameter pair ((In ^-,0), 0). 
When we choose a long enough existence interval [0, s R ] such that the I value at 
s R of the above solution is less than ? — 1, we have K L {oS) < <r — 1 as u) is close 
to enough by lim l* R ito) = In above and the "continuous dependence" 

property. This contradicts to lim 1* l {oj) = Thus lim 1* l {lo) = — oo and 

then (iii) follows. ■ 
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Note that by Lemma El and its proof we can roughly say that given u)\ < 
lu 2 , (I(s;uj2),v(s;uj2)) traces out an orbit curve which is closer to (ln-^-,0) 
than (1(s;lui),v(s;lu 1 )) does as s — > oo, i.e., (l(s; 0J2), v(s; LU2)) is "enclosed" 
by (l{s)U X ),v{s\u)x)). 

Lemma 7 Given any 1(0) E (-00, In 2gi), G, andie {0,1,2,...}, th ere is at 

most one to = u>(G, i,l(0)) among (0, ^f) such that f?3 p ,(T?P, f75J) ftas a 5 /ofeaZ 
solution (1(-;lu),v(-;lu)) on [0, 00) which has exact 2i intersection points with the 
I— axis as s increases from to 00. 

Proof. Note that any global solution (1(-;lo),v(-; u>)) of 113(1 .1(14 (1 has the same 
orbit curve as (1(-;lu),v(-;lu)), i.e., (l(s; lu),v(s;lu)) = (l(s + sq;lu),v(s + So;u>)) 
for some so by Proposition [5j Therefore without loss of generality, we consider 
(/(•; lu), «(•; ui)) on [0, 00). Fix I = 1(0) and G in 1(15(1 . Note that the initial value 
(1(0), v(0)) is on the curve: 

v = I(1;lu,G) : =-^ e -< + ^. (26) 

On one hand vq := v(0) = 1(1(0); ■, G) is a decreasing function of lu on (0, 00). 
Let (li(-;w), Vi(-; to)) be the segment of (l(-;ui),v(-;u)))'s orbit curve connect- 
ing Iu+i)ji(w) and 1*l(oS) in the lower I — v plane. Then on the other hand, 

given u>\ < u>2 and S\,S2 such that li(si;u){) = h(s2]U2) = '(0), it can be eas- 
ily seen from Lemma (i) and its proof that vn(s\\u)i) < Vi(s2',LU2), if both 
Vi(si;u>i),Vi(s2;iJ2) are nonpositive, i.e., increasing in lu. Since there is at most 
one intersection point for increasing and decreasing functions, this lemma then 
follows. ■ 

Now we state our main result: 

Theorem 8 (i). Given any 1(0) 6 (— oo,ln2gi) and i 6 {0,1,2,...}, there is 
a corresponding Gi = Gj(/(0)) > such that for any given G € [0, Gj] 

there is a unique lu = lu (G;i,l(0)) £ (0,^-) such that ff3|) . {77|). 
/ias a global existence solution (1(-;lo),v(-;lj)) on (—00,00) which rotates 
clockwise around (ln-^-,0), /ias 2i or 2i + 1 (exact 2i for G — case) 
intersection points with the I— axis as s increases from to 00, and after the 
2i or2i + l intersections we have v(s; to) < 0, l(s; lu) \ —00, v(s; lu) — > 0~ 
as s > 00. Moreover, given any G € fl[0, Gi(l(0))), lu(G; -,l(0)) is a 

i 

decreasing function with the properties that Lu(G;i,l(0)) \ + as i — > 00 
for 1(0) e (-00, In 2gi)\{ln G <E n[0, G,(Z(0))) and w(0; i, In \ 0+ 

as i — ► 00. 

(«i). For i(0) G [In 2^,00) and any G, we Ziaue that lu = is the only num- 
ber such that \13ti , fT^P , if J5|) /ias global existence solution (l(s; 0), v(s; 0)). 
Moreover, if G > 0, then there is some d S (0, 00) swc/i i/iai u(s; 0) > 
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for s £ [0, •&), 0) = 0, and v(s; 0) < for s £ oo) with the asymp- 
totic behavior /(s;0) \ — oo, v(s;0) — ► 0~ as {) < s — ► oo; if G < 0, 
i/ien w(s;0) < /or s G (0, oo) iwi/i i(s;0) \ — oo, v(s;0) — > 0~ as 

< s — ► oo. 

(Hi). If G £ (— oo,— Vq], IJ.'j|) , f77| ), i|J.5|) /ias no global existence solutions. When 
G £ {— Vo,0), given any 1(0) £ (m[ ^-jf J,m[ J) ana 

1 £ {0,1,2,...}, there is a unique uj — u(i,l(0);G) £ (0, -jj-) sucft i/iai 
y3^^2P,y3j) ^ as a 9^°^ existence solution (1(-;lo),v(-;uj)) satisfying all 
the properties described in (i) ( but having exact 2i intersection points 
with the I— axis as s increases from to oo here); on the other hand, 

i)-?5J) ,i fm ), i)-?5|) has no global existence solutions whenever 1(0) (ln[— — 
in this case. 

First, we explain the physical meaning of the above results and then prove 
this theorem. 



3.1 Interpretation of Theorem |S] 

(i). Remember that 1(0) = \jxkq, where no is the tip's curvature and G is 
the tangential velocity of the tip. Theorem [S] (i) tells us that given a weakly 
excitable medium (V$,D are fixed), if the tip's curvature «o is not too large, 
i.e., Ko G (0, ^ )j there is a corresponding range for the tangential velocity 
of the tip so that for each G in this range there is a unique rotating frequency 

co(G; i, hiKo) £ (0, -p 2 -) such that this medium supports a normally propagating, 
steadily rotating, and growing (since G > 0) plane curve with angular frequency 
ui(G; i, ln^o) and tip curvature kq>, the curvature of which n(s) — e'w changes 
its monotonicity 2i or 2i + 1 times (exact 2i for G — case) as arc length s — > oo 
and finally decreases to in the Archimedean spiral's way, i.e., lim k(s) 2 s = 

s — >oo 

u(G;iM K0 ) ^ gee ^ . n j2|j^ Especially; the case «j = 0, G = 0" is just the one 
appearing in |75) . 

(ii) . On the other hand, if the tip's curvature kq is too large, i.e., kq £ 
[2^-,oq), then the medium supports only ll nonrotating^ spiral waves (to = 0). 
Moreover, given any G > (G < 0, respectively), there is a unique growing 
(contracting, respectively) spiral wave with tip's curvature n and tangential 
velocity G, the curvature of which changes its monotonicity exactly once (0 
time, respectively ) and the asymptotic behavior of which is the same as in 
above (i) but the curvature k(s) decays to in different order, not Archimedean 
one (see THEOREM 2 of [22). 

(iii) tells us that it is impossible for an excitable medium to support a spiral 
wave with too large contracting velocity of the tip, i.e., G £ (— oo,— Vq]. For 
smaller contracting velocity, G £ (— Vb,0), there is a corresponding curvature's 

range of the tip, ( — ^3 G - , ^ 0+ aAo Z ^LLy such that when the tip's curva- 
ture kq does not belong to this range, there is no any spiral wave but when given 



CLASSIFICATION OF STEADILY ROTATING SPIRAL WAVES 



15 



any kq in this range, there is a unique u)(i, In kq; G) € (0, — jj-) such that this 
medium supports a spiral wave with tip's contracting velocity G, tip's curvature 
kq. rotating frequency ui(i, In kq; G), and the same properties as mentioned in 
(i) above (but the curvature changes its monotonicity exact 2i times). 

Note that we also have similar explanations as above by replacing k, u>, Vq 
with — k, —u), — Vq, respectively when k < 0. 



3.2 Proof of Theorem |U 

Proof. As explained in the proof of Lemma0 it suffices to consider (/(•; cu), «(•; cj)), 
Yl 

D 



Letting uj = (JzoJ becomes 



^/(^,G):=-2(^) 2 e- ; + §. (27) 
Plugging (|27[) into E(l,v), we obtain 

S(l,-2$)»e-' + ||) ^ ^ 2i [4(^) 4 + e 4i -^ + (§) 2 e 2 '-4(§)(^)V]. 

(28) 

It can be easily derived that lL x A — 2ax i + 4a 4 > 0,Vx £ R,a j£ 0" which implies 
£(Z,-2(^) 2 e-') > for all I e R (i.e., when G = 0). Hence given any I, we 

have that /(Z; -jj-,0) <the negative t> value of curve u E(l,v) ~ 0", i.e., curve 
lf2*7jl lies below E(l, v) = when G = 0. Then by continuity, there is a G(Z) > 

such that 1(1; ^jj-, G) < the negative v value of curve "E(l, v) = 0", whenever 
G G (— oo, G(l)]. From above observation and Lemma|3| we have that if there is 

some lo > such that (fT3|) . (fl~H) . fl5|) has global solutions, then uj G (0, ^jf-)- Thus 
we have obtained the uniqueness property by Lemma[7| It is left for (i) to prove 
the "existence" of such uj. The strategy is to find suitable uj such that there are 
intersections points between the curve l|26(l and some special segment of the orbit 
curve of (1(-;uj),v(-;uj)), which implies (|13|I . (|14(I . H15|I has global solutions with 
the properties mentioned in this theorem. By Lemma[Sl the ranges of ^(j +1 )p(') 

and l* L (-) over (0,^-) are open intervals, say, (iZj+i, In 2^-) and (— oo, Li) 
respectively, where Iq L (-) '■= — oo and Lq :— — oo. For any i G {0,1,2,...} 

and 1(0) G (Li,R i+1 ), Vi(l(0);uj) is well-defined for any lo G (0,^-), where 

Vi(l;uj) := v(s;ui) if there is some s such that 1(s;uj) = I and (l(s; uj), v(s; uj)) 
is on the "lower" orbit curve connecting ^* i+1 )p(w) and 1* l (uj). Thus given 
any 1(0) G (L. t ,R i+1 ) and G G (-oo, G;(Z(0))], where Gj(-) := G(-), there is 

some positive close to such that Vi(l(0); cjpj >the negative u value of 
"£(Z(0),v) = 0" > I(l(0);u) Pi ,G) by Lemma the continuity of 1(1(0); G), 
and the observation in the beginning of this proof; for 1(0) = Li or R4+1, it 
still holds by Lemma El - Now consider 1(0) G (i? i+ i,ln^) (or 1(0) G 

(— 00, Li)). By definition, there is uj G (0, ^-) such that l( i+ i) R (^o) — '(0) (or 
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l* L (u>o) = 1(0))- We require now that 1(1(0); uiq, G) < which is equivalent to 

G<tu e- m . (29) 

Let Gi(l) := oJoe~ l . Therefore, given any G £ (— oo, Gj(Z(0))], there is also some 
uop i close to ujq such that Vi(l(0); wpj > 1(1(0); uip t , G) by virture of Lemma 
(i),(ii). On the other hand, given any G S [0, oo) and 1(0) G (— oo, ln^ 1 ), 
there is some positive close to such that Ut(i(Q);wjv 4 ) is well-defined 
and Wi(?(0);wAr.) < 7(/(0); G), since lim 7(1(0); w,G) = £ > and 

lim Vi(l(0);uj) =the negative w value of u E(l(0),v) = 0" < by Lemma El 

(iii) and the continuous dependence theorem of o.d.e. Note that Vi(l(0);-) is 
a continuous function on [wjv^wpj by Lemma |H] (ii) . Obviously 7(Z(0);-,G) is 
also a continuous function. Therefore choosing 1(0) £ (— oo,ln-^) and then 

G G [0, Gi(/(0))), there must be some <7j € (wjv^wpj C (0, ^jy-) such that 
Ui(/(0);(Ti) = I(l(0);ai,G) by the intermediate value theorem of continuous 
functions. Let ui(G;i,l(0)) := cr^. Hence "existence" is guaranteed. Given 
G € n[0, G,(i(0))), the decreasingness of cu(G; -,1(0)) as i — > oo immediately 

follows from the explanation after Lemma|Hl i.e., u (l(s;u>2),v(s;u>2)) is enclosed 
by (l(s; uii), v(s; u>i)) as lo\ < uj 2 " and the decreasingness of 7(2(0); •, 0). Now we 
prove ]imu(G;i,l(0)) = for 1(0) £ (-oo, In ^)\{ln 3.}, G £ n[0, Gj(/(0))) 

and lim w(0; z, In ^) = 0. For the first part, assume the contrary that lim uj(G; i, 1(0)) = 

u f (G;l(0)) > 0. Note that lim l* R (u>f(G; 1(0))) = lim l* L (tu f (G; 1(0))) = In 3 

i — >oo i — >oc 

By 1(0) ^ In there must be some i large enough such that 1(0) £ (l^ i+1 - )R (ujf(G;l(0))),ln 

C ( J R i+1 ,ln^)orZ(0) £ (-oo, l* L (ujf(G; 1(0)))) cj-oo, L,). Arguing as before, 
then there is some to £ (0,Lu f (G;l(0))) such that Vi(l(0);u>) = I(l(0);u>,G). By 
uniqueness, this to must be cu(G;i,l(0)) which leads to a contradiction, since 
it is impossible that uj(G;i,l(0)) < ujf(G;l(0)). For the second part, note that 
Ui(ln-^;a;/(0;ln-^-)) is always defined and lim Vi(ln ^; ujf(0; In jj)) = as ex- 

plained in the last of the proof of Lemma Thus for i large enough we have 
^(ln^; W/ (0;lnf )) > /(ln$; W/ (0;m§ ),0), since /(In f ; to f (0; In f ), 0) < 
0. Then again by Lemma[S](iii) and the intermediate value theorem, there must 
be some lj £ (0, ujf(0; In jy)) such that Vi(ln^-;uj) = /(In ^; u, 0) which also 
leads to a contradiction. This completes (i). 

Part (ii) follows from LemmaEDand THEOREM 2 in [23J immediately. 

Now we prove (iii). When G £ (— oo, — Vq], the curve (I, I(l;co,G)) entirely 
lies outside the curve "E(l, v) = 0" for any u> > 0. Thus by Lemma|3 we have 
that lfT3|) . ((T11l . ((r5j) has no global existence solutions. When G £ (-Vo,0), the 

parts of the curve (/, 1(1; w,G)), where 1(0) £ (ln[ (y °~^f ~ G -^], ln[ (y ° + ^ ^]), 
also lie outside the curve u E(l,v) — 0" for any uo > 0. As explained above, we 

only need to consider 1(0) £ (ln[ (Vb ~^f ~^ ], m [iMvjpEg!)]) now . Note 
that the right-hand side of (|28H is greater than V I, for any G < and (|29|l is 
also valid. Arguing as in (i), we then have completed (iii). 
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Therefore we have finished the proof. ■ 



4 Numerical Results 

The following formulas can be referred to |4(J1 p. 45]. 
Let 

O (t) =wi + 0o(O) (30) 

be the angle of the tip at time t, where u> is the rotating frequency determined 
by Theorem [8] and 0q(0) is the initial angle of the tip, which can be chosen 
arbitrary. Let also 9(s,t) be the angle of the position away from the tip at a 
distance "arclength s" at time t. Then 

9(s,t) = 6 (t) - [ e l ^d P , (31) 



where l(s) is the function from Theorem [S] 

The cartesian coordinates of the tip at time t are given by 

-^X (t) = (De m - V ) sin 6 {t)-G cos 6 {t) = (De m - V ) sin[orf + O (O)] - Gcos[u;t + 9 (0)} and 

(32) 

4^Y (t) = (Vq - De m ^j cos9 {t) - Gsm6 a (t) = (v - De m ^j cos[wt + O (O)] - G sm[ut + 6 (0)], 

where the initial condition (Xq(0), Yq(0)) can be chosen arbitrarily. Then, the 
cartesian coordinates (X(s, t), Y(s, t)) of the position away from the tip at a 
distance "arclength s" at time t are 

X(s,t)= [ cos9(Z,t)d£ + X (t) = [ cos[ujt + 9 (0) - e l( - p} dp]d£ + X {t) 
Jo Jo Jo 

(33) 

Y(s,t)= I sm0(Z,t)d4 + Y o (t)= [ sin[wt + o (O)- / e l ^dp]d$ + Y (t), 
Jo Jo Jo 

which is the same as (4.43) in p. 45]. 

Now we explain how to determine the parameters D, Vq, 0o(O), 1(0), G, w 
and then to draw the pictures of spiral waves in Theorem |H1 

First, parameters D > 0, Vq > 0, 0o(O) can be chosen arbitrarily. 

In Theorem El (i), we choose 1(0) G (-oo,ln^), i £ {0,1,2,3,...} first. 
Then we choose G G [0, G,(Z(0))] (for example choose G = 0). Finally, there is 

a unique u — uj(G; i, 1(0)) £ (0, ^jy-) by Theorem (i). Plug above parameters 
into il-il'l) and (|33|l and then the corresponding spiral wave has the properties 
described in Theorem (i). 

For (ii), we choose 1(0) £ [In =ff, oo), w = 0, and any G. Plug above parame- 
ters into (|32|l and l|33l) and then the corresponding spiral wave has the properties 
described in Theorem |H1 (ii). 
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For (iii), we first choose G € {-V , 0) and then 1(0) e (ki[ — ^ ° !2],ln[lM^zg!2]) 

and any i £ {0, 1, 2, 3, ...}. Finally there is a unique ui(i, Z(0); G) G (0, ^jj-)- Also, 
plugging above parameters into l|32(l and (|33|l . the corresponding spiral wave 
then has the properties described in Theorem |S| (iii). 

Acknowledgements C.-P.Lo thanks Mr. J.-C.Tsai and professors J.-S.Guo 
and S. -S.Lin for helpful suggestions. 
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